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I. INTRODUCTION 

A significant range of studies of ultracold atomic sys- 
tems have focused on their complex quantum behavior 
in various interaction processes. Among these, special 
attention has been payed to the quantum interface be- 
tween light and matter, and quantum memory in partic- 
ular PHl . Most of the schemes for light storage in atomic 
ensembles are based on idea of the A-type conversion of 
a signal pulse into the long-lived spin coherence of the 
atomic ground state. The electromagnetically induced 
transparency (EIT) protocol in a warm atomic ensem- 
ble was successfully demonstrated in Ref. and also 
in Ref. (EJ, where a single photon entangled state was 
stored in two ensembles of cold atoms with an efficiency 
of 17%. Recent experiments on conversion of a spin po- 
lariton mode into a cavity mode with efficiency close to 
90% @ and on the narrow-bandwidth biphoton prepara- 
tion in a double A-system under EIT conditions Q show 
promising potential for developing a quantum interface 
between light and atomic systems. However, further im- 
provement of atomic memory efficiencies is a challenging 
and not straightforward experimental task. In the case of 
warm atomic vapors, any increase of the sample optical 
depth meets a serious barrier for the EIT effect because 
of the rather complicated and mainly negative rule of 
atomic motion and Doppler broadening, which manifest 
in destructive interference among the different hyperfinc 
transitions of alkali- metal atoms |8| . In the case of ultra- 
cold and dilute atomic gas, which can be prepared in a 
magneto-optical trap (MOT), for some experimental de- 
signs optical depths around hundreds are feasible but 
there are a certain challenges in accumulating so many 
atoms and making such a system controllable. One pos- 
sible solution requires special arrangements for effective 
light storage in MOT in a diffusion regime; see Ref. [Toj . 

Recent progress in experimental studies of light local- 
ization phenomenon in the dense and strongly disordered 



atomic systems (111 . 1 1 21 ] encourages us to think that the 
storage protocols for light could be organized more ef- 
fectively if atoms interacted with the field cooperatively 
in the dense configuration. If an atomic cloud contains 
more than one atom in the volume scaled by the radi- 
ation wavelength, the essential optical thickness can be 
attained for a smaller number of atoms than it is typically 
needed in dilute configuration. In the present paper we 
address the problem of light scattering by such an atomic 
system, which has intrinsically cooperative behavior. Al- 
though the problem of cooperative or dependent light 
scattering and super-radiance phenomenon have been 
well established in atomic physics and quantum optics 
for decades (see Refs. [H, [3), microscopic analysis for 
the atoms with degenerate ground states is still quite 
poorly performed in the literature floT |. The microscopic 
calculations reported so far have been done mostly for 
"two- level" atoms and were basically motivated by the 
problem of mesoscopic description of the light transport 
through disordered media and by an Anderson-type lo- 
calization, where transition from weak to strong disorder 
plays a crucial role; see Refs. [l6| - [l8j . 

In this paper we develop a microscopic theory of 
the cooperative light scattering from an atomic sys- 
tem consisted of A-type configured atoms with the spin- 
degenerate ground state. The results are compared with 
a similar system of "two-level" atoms of the Dicke model. 
We discuss advantages of the considered system in the 
context of its possible implications for the problem of 
light storage in a macroscopic ensemble of dense and ul- 
tracold atoms. 



II. THEORETICAL FRAMEWORK 

A. Transition amplitude and the scattering cross 
section 
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The quantum-posed description of the photon scatter- 
ing problem is based on the formalism of T matrix, which 
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is denned by 

f(E) = V + V — — ^V, (2.1) 
E-H 

where H is the total Hamiltonian consisting of the non- 
perturbed part Hq and an interaction term V such that 
H = Hq + V. The energy argument E is an arbitrary 
complex parameter in Eq. (|2.1|) . Then the scattering pro- 
cess, evolving from initial state \i) to the final state |/), is 
expressed by the following relation between the differen- 
tial cross section and the transition amplitude, given by 
the relevant T-matrix element considered as a function 
of the initial energy Ec 

dVi^f = 7?77^2 l^e'k', 9 ek(^ + i0)| 2 dn (2.2) 
fi c 4 (27Tj z 

Here the initial state \i) is specified by the incoming pho- 
ton's wave vector k, frequency u> = lo^ = ck, and po- 
larization vector e, and the atomic system populates a 
particular ground state \g). The final state |/) is speci- 
fied by a similar set of the quantum numbers, which are 
additionally upscribed by the prime sign, and the solid 
angle f2 is directed along the wavevector of the outgoing 
photon k'. The presence of quantization volume V in this 
expression is caused by the second quantized structure of 
the interaction operators; see below. The scattering pro- 
cess conserves the energy of input and output channels, 
such that Ei — Ej. 

Our description of interaction process of the electro- 
magnetic field with an atomic system is performed in 
the dipole approximation. This states that the origi- 
nal Hamiltonian, introduced in the Coulomb gauge and 
valid for any neutral charge system, has been unitarily 
transformed to the dipole-type interaction with the as- 
sumption that atomic size is much smaller than a typical 
wavelength of the field modes actually contributing to 
the interaction dynamics. Such a long-wavelength dipole 
approximation see Ref. for derivation details leads to 
the following interaction Hamiltonian for an atomic en- 
semble consisting of N dipole-type scatterers interacting 
with the quantized electromagnetic field: 

N 

V = -^dWE(r )+4eif, 

a=l 

N 

= EyEH w ) (2-3) 

a—1 s 

The first and most important term is normally inter- 
preted as interaction of an ath atomic dipole with 
electric field E(r) at the point of dipole location. How- 
ever, strictly defined in the dipole gauge, the latter quan- 
tity performs the microscopic displacement field, which 
can be expressed by a standard expansion in the basis 
of plane waves s = k, a (where a — 1,2 numerates two 
orthogonal transverse polarization vectors e s = ek a for 



each k) 

E(r) = E( + '(r)+E(-)( r ) 

= J2 (^) V2 [ie s a s ^ - ie.ate"*"] 

N 

= E ± (r) + Y E e s (e s dW)e*M r -^ (2.4) 

6=1 s 

Here a s and a\ are the annihilation and creation opera- 
tors for the sth field's mode and the quantization scheme 
includes the periodic boundary conditions in the quanti- 
zation volume V. The bottom line in Eq. (|2.4[) indicates 
the important difference between the actual transverse 
electric field denoted as Ej_(r) and the displacement field. 
The difference cannot be ignored at the distances compa- 
rable with either atomic size or the radiation wavelength, 
which is the subject of the present report. For such a 
dense configuration the definitions (|2.3[) and (12. 4[) should 
be clearly understood. 

Let us make a few remarks. The second term in 
Eq. (I2.3[) reveals a nonconverging self-energy (self-action) 
of the dipoles. This term is often omitted in practi- 
cal calculations since it does not principally affect the 
dipoles' dynamics, particularly when the difference be- 
tween transverse electric and displacement fields is small. 
It can be also formally incorporated into the internal 
Hamiltonian associated with the atomic dipoles. How- 
ever, as was pointed out in Ref. [l8| via tracing the 
Heisenberg dynamics of atomic variables, the self-action 
term is mostly compensated by the self-contact dipole 
interaction. The latter manifests itself in the dipoles' dy- 
namics when r = r a = rh for interaction of a specific 
a-th dipole in Hamiltonian (|2.3[) with the longitudinal 
field created by the same dipole in the second term in 
Eq. (|2.4p . Both these nonconverging self-action and self- 
contact interaction terms can be safely renormalized in 
evaluation of a single-particle contribution into the self- 
energy part of the perturbation theory expansion for the 
resolvent operator; see below. 



B. Resolvent operator and TV- particle Green's 
function 

The transition amplitude (I2.1[) can be simplified if we 
substitute in it the interaction operator (|2.3p by keeping 
only the terms with annihilation of the incoming photon 
in the input state and creating the outgoing photon in the 
output state. Such a simplification is in accordance with 
the standard approach of the rotating wave approxima- 
tion, which is surely fulfilled for a near-resonance scatter- 
ing process. As a consequence of this approximation the 
transition amplitude is now determined by the complete 
resolvent operator projected onto the vacuum state for 
the field subsystem and onto the singly excited state for 
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the atomic subsystem 

R(E) = PR(E)P = P- 



n, r a and b' 



1 



-P. 



E-H 

Here we denned the following projector 

N 

|mi, . . . ,m a -i, n,m a +i, ■ ■ 



(2.5) 



E E 



■ m N ) 



^ = E 

o=l {ri 

(mi,...,m Q _i,n, m Q+ i,...,mjv| x |0)(0| Field (2.6) 

which selects in the atomic Hilbert subspace the entire 
set of the states where any jth of N — 1 atoms populates 
a Zeeman sublevel \rrij} in its ground state and one spe- 
cific oth atom (with a running from 1 to N and j a) 
populates a Zeeman sublevel \n) of its excited state. The 
field subspace is projected onto its vacuum state and the 

operator R(E) can be further considered as a matrix op- 
erator acting only in atomic subspace. The elements of 
the T matrix can be directly expressed by the resolvent 
operator as follows: 



T s Vk', S ek(£) 



V 



E 



b,a—l n'n 



(de'):, K (de)„ ma e-*' r ^ kr « 
.m' b _ 1 ,n',m' b+1 ...\R(E)\ 



. .m a -i,n, m a+ i . . .) 

(2.7) 



This performs a generalization of the well-known 
Kramers-Heisenberg formula [2(| for the scattering of a 
photon by a many-particle system consisting of atomic 
dipoles. The selected specific matrix element runs all 
the possibilities when the incoming photon is absorbed 
by any ath atom and the outgoing photon is emitted by 
any &th atom of the ensemble, including the possible co- 
incidence a = b. The initial atomic state is given by 
\g) = | mi, . . . , TOzv) and the final atomic state is given by 
W) = \m' 1 ,...,m' N ). 

The projected resolvent operator contributing to Eq. 
()2.7p is defined in the the Hilbert subspace of a finite size 
with dimension d e N d^ -1 , where d e is the degeneracy of 
the atomic excited state and d g is the degeneracy of its 

ground state. The matrix elements of operator R(E) can 
be linked with the iV-particle causal Green's function of 
atomic subsystem via the following Laplace-type integral 
transformation: 

(. . .m'^^n' ,m' b+l . . . \R{E)\ . . .m a -i,n,m a+ i . . .) 

x6 (ri - rx) . . . S (r' 6 - r b ) . . . S (r' tt - r tt ) . . . 5 (r' N - r N ) 



ft Jo 



dt exp 



-Et 



G^(l',t;...;b',t;...;N',t\l,0;...;a,t;...;N,0)(2.8) 

where on the right side we denoted j = rrij , Tj (for j =/= a) 
and j' — TripTj (for j' ^ &), and for specific atoms a = 



n' , v' b . Here = r^-, for any j = 1 -j- N, is 
the spatial location of jth atom, which is assumed to be 
conserved in the scattering process. This circumstance 
is expressed by the sequence of /delta functions in Eq. 
(|2.8p . The causal Green's function is given by the vacuum 
expectation value of the following chronologically (T)- 
ordered product of atomic second quantized VP operators 
introduced in the Heisenberg representation 

G^(l\t' 1 ;...;b\t' b ;...;N\t' N \l,t 1 ;...;a,t a ;...;N,t N ) 
= (T\f m; (ri ,t[)... tt n , (r' 6) t' b ) ...9 m , N (r' N , t' N ) 

( r w, t N ) . . . *t (r a , t a )... (ri.ii)), 

(2.9) 

where ^...(. . .) and (...) are respectively the annihi- 
lation and creation operators for an atom in a particu- 
lar state and coordinate. All the creation operators in 
this product contribute to transform (|2.8[) while being 
considered at initial time "0" and all the annihilation 
operators are considered at a later time t > 0. That al- 
lows us to ignore effects of cither bosonic or fermionic 
quantum statistics associated with atomic subsystem as 
far as we neglect any possible overlap in atomic loca- 
tions and consider the atomic dipoles as classical objects 
randomly distributed in space. We ordered operators in 
Eq. (I2.9[) in such a way that in the fermionic case (un- 
der the anticommutation rule) and without interaction it 
generates the product of independent individual single- 
particle Green's functions associated with each atom and 
with positive overall sign. 

The perturbation theory expansion of the iV-particlc 
Green's function (|2.9[) can be visualized by the series of 
the diagrams in accordance with the standard rules of 
the vacuum diagram technique; see Ref. 20]. After re- 
arrangement the diagram expansion can be transformed 
to the following generalized Dyson equation: 




b=a 
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where the long straight lines with arrows correspond with 
individual causal single-particle Green's functions of each 
atom in the ensemble such that the first term on the 
right side performs the graph image of nondisturbed N- 
particle propagator (|2.9p . The dashed block edged by 
short lines with arrows is the complete collective TV- 
particle Green's function dressed by the interaction. In 
each diagram block of equation (|2.10j) we indicated by 
a, 6, c (running from 1 to N) the presence of one specific 
input as well as an output line associated with the single 
excited state equally shared by all the atoms of the en- 
semble. The sum of tight diagrams, which cannot be re- 
duced to the product of lower order contributions linked 
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by nondisturbed atomic propagators, builds a block of so- 
called self-energy part E. The diagram equation (12.10)) 
in its analytical form performs an integral equation for 
G( N \. . .). With its transformation to the energy repre- 
sentation (|2.8[) the integral equation can be recomposed 
to the set of algebraic equations for the matrix of the 

projected resolvent operator R(E), which can be further 
numerically solved. The crucial requirement for this is 
the knowledge of the self-energy part (quasi-energy op- 
erator acting in the atomic subspace), which as we show 
below can be approximated by the lower orders in expan- 
sion of the perturbation theory. 



should be associated with the longitudinal self-contact 
interaction. These terms are compensated by the dipo- 
lar self-action; see Eq. (|2.3p and the related remark given 
above. The residual nonconvergency has radiative nature 
and demonstrates general incorrectness of the Lamb-shift 
calculation in assumptions of the long-wavelength dipole 
approximation. Finally we follow the standard renormal- 
ization rule, 



^l(E) = ^(E)S n , n , 



ih- 



(2.14) 



C. The self-energy part 

In the lower order of perturbation theory the self- 
energy part consists of two contributions having single- 
particle and double-particle structures. Each specific line 
in the graph equation (|2.10p associated with excitation 
of an ath atom generates the following irreducible self- 
energy diagram: 




where Al — > oo is incorporated into the physical energy 
of the atomic state. To introduce the single-atom nat- 
ural decay rate 7 we applied the Wigner-Weiskopf pole 
approximation and substituted the energy E — Tilu^. + E g 
by its near resonance mean estimate E w E n with as- 
sumption that the atomic ground state is the zero-energy 
level such that E g = Ylj=i E mj = E m = 0. Then energy 
of the excited state is given by E n = hujQ, where ujq is 
the transition frequency. 

In the lower order of perturbation theory, the double- 
particle contribution to the self-energy part consists of 
two complementary diagrams: 



(2.11) 



1 



E — hu> — E„ 



iO 



which is analytically decoded with applying transforma- 
tion (12.81) in the energy representation. Here the inter- 
nal wavy line expresses the causal-type vacuum Green's 
function of the chronologically ordered polarization com- 
ponents of the field operators 



D$ (R,r) = H(r' ) i')i(r,^ 1 



(2.12) 



which depends only on difference of its arguments R = 
r' — r and r = t' — t and has the following Fourier image: 

/oo 
-OO 



U\ 3 f.2, (1) / a.' 



3° \^ R ) 6 ^ 



X^Xy 1 

~ 3 ^ 



see Ref. |20|. Here h£\...) with L = 0, 2 are the spheri- 
cal Hankel functions of the first kind. As follows from Eq. 
(|2.11j) the Green's function (|2.13l) contributes in that ex- 
pression in a self-interacting form with spatial argument 
R —> 0. As a consequence the expression (|2.11[) becomes 
non-converging in the limit R — > and the integration 
over uj is nonconverging. Part of nonconverging terms 



(2.13) 




did 

2^ 



1 



E — Hoj — E m - E, m i + iO 



( 2 - 15 ) 



and 




1 



E + huj-E n - E n > + iO 



which are responsible for the excitation transfer from 
atom a to atom b separated by a distance R a b- The vec- 
tor components of the dipole matrix elements d v m , n and 
dn'm are related with atoms a and b respectively. In the 
pole approximation E ss E n = Kuiq the 6 function fea- 
tures dominate in the spectral integrals (I2.15[) and (|2.16[) 
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and the sum of both the terms gives 

K mV;nm( £ ) ~ S roV;nm (^o) + ^m^n']nm (^0 ) 

= ^<, m C'„4^( R ^^o). (2.17) 

The derived expression has clear physical meaning. For 
nearly located atoms the real component of the double- 
particle contribution to the self-energy part reproduces 
the static interaction between two atomic dipoles. Its 
imaginary component is responsible for formation of co- 
operative dynamics of the excitation decay in the entire 
radiation process. For long distances, when the atomic 
dipoles are separated by the radiation zone, this term 
describes radiation interference between any pair of two 
distant atoms, which weakly reduces with the interatomic 
separation. For short distances or in a dense sample the 
cooperative effects become extremely important and the 
scattering process becomes strongly dependent on a par- 
ticular atomic configuration. 

It is a challenging problem to further improve the self- 
energy part by taking into consideration the higher or- 
ders of the perturbation theory expansion. Here we only 
substantiate the validity and sufficiency of the lower or- 
der approximation for the considered configuration. The 
main physical reason for this is weakness of interaction. 
This justifies ignoring of any deviation from free dynam- 
ics of atomic variables on a short time scale associated 
with the light retardation on distances of a few wave- 
lengths. That yields main cooperation in the radiative 
dynamics among neighboring dipoles which can effec- 
tively interact via static longitudinal electric field. The 
diagram (|2.16p . in contrast with (|2. 15|) . is mostly im- 
portant for evaluation of the static interaction such that 
in this graph the field propagator preferably links the 
points with coincident times on atomic lines. As a con- 
sequence, the presence of such diagram fragments as a 
part of any irreducible diagrams in higher orders would 
make the overall contribution small and negligible just 
because the static dipolc-dipole interaction only weakly 
affects the dipoles' dynamics during the short retardation 
time, which can be roughly estimated by the wave period 
27r/wo- For the same reason we can ignore any vertex- 
type corrections to the diagram (|2.14|) . Another part of 
the self-energy diagrams in higher orders can be associ- 
ated with correction of the static interaction for itself. If 
the atomic system were as dense as the atoms were sep- 
arated by a distance comparable with atomic size (much 
shorter than the radiation wavelength) then the descrip- 
tion of the static interaction in the simplest dipole model 
would be inconsistent and insufficient. This correction is 
evidently ignorable for atomic ensemble with a density 
of a few atoms in a volume scaled by the cubic radiation 
wavelength. In this case the higher order static correc- 
tions are negligible as far as the dipole-dipole interaction 
is essentially less than the internal transition energy. As 
we can finally see, for the considered atomic systems, the 
self-energy part is correctly reproducible by the intro- 
duced lower order contributions. 
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FIG. 1: (Color online) The excitation diagram of "two-level" 
V-type atom (left) and overturned "two- level" A- type atom 
(right). In both the configurations the light scattering is con- 
sidered for the left-handed <j_ polarization mode. The A-atom 
populates the Zeeman sublevel with the highest angular mo- 
mentum projection. 



III. RESULTS AND DISCUSSION 

A. Cooperative scattering from the system of two 
atoms 

Let us apply the developed theory to the calculation 
of the total cross section for the process of light scatter- 
ing from the system consisting of two atoms. We con- 
sider two complementary examples where the scattering 
atoms have different but similar Zeeman state structure. 
In the first example we consider V-type atoms, which 
have Fo = total angular momentum in the ground state 
and F = 1 total angular momentum in the excited state. 
Such atoms are the standard objects for discussion of 
the Dicke problem see Ref. [l3|], and each atom per- 
forms a " two- level" energy system sensitive to the vector 
properties of light. In an alternative example we con- 
sider the A-type atoms, which can be also understood as 
overturned "two- level" system, which have Fo = 1 to- 
tal angular momentum in the ground state and F = 
total angular momentum in the excited state. For the 
latter example in the scattering scenario we assume the 
initial population by atoms of a particular Zeeman sub- 
level of the ground state, which has highest projection 
of the angular momentum. Both the excitation schemes 
and transition diagrams in the laboratory reference frame 
are displayed in Fig. [TJ 

In Figs. [2E1 we reproduce the spectral dependencies 
of the total cross section for a photon scattering from 
the system consisted of two atoms separated by differ- 
ent distances R and for different spatial orientations. 
The variation of interatomic separation from R = 10A 
(independent scatterers) to R — 0.5A (strongly depen- 
dent scatterers) transforms the scattering process from 
its independent to cooperative dynamics. In the plotted 
graphs the frequency spectra, reproduced as function of 
the frequency detuning A = lo — luq of the probe fre- 
quency uj from the nondisturbed atomic resonance wo, 
are scaled by the natural radiation decay rate of a sin- 
gle atom 7, which is significantly different for A- and 
V-type energy configurations, such that 7(A) = ^(V). 
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As a consequence the near-field effects responsible for 
the resonance structure of the resolvent operator and the 
cross section manifest more perceptibly for the V-type 
atoms, which are traditionally considered in many dis- 
cussions of the Dicke system in literature. In the sym- 
metric collinear excitation geometry, when the internal 
reference frame coincides with the laboratory frame see 
Fig. [5] the left-handed <r_ excitation channel shown in 
Fig. [T]is only allowed for either V- or A-type transition 
schemes. In such a symmetric configuration the inter- 
atomic interaction via the longitudinal as well as via the 
radiative transverse fields splits the excitation transition 
in two resonance lines. For the case of the V-type excita- 
tion the observed resonances demonstrate either superra- 
diant or subradiant nature. This is an evident indicator 
of the well-known Dicke effect of either cooperative or 
anticooperative contribution of the atomic dipoles into 
the entire radiation and scattering processes; see Refs. 
[l3l [l9j . For the system of two A-type atoms separated 
by the same distances the atomic line also splits in two 
resonances, but they are less resolved and have relatively 
comparable line widths. The spectral widths indicate 
a slight cooperative modification, which is much weaker 
effect than in the case of V-type atoms. The physical rea- 
son of that is the contribution of the Raman scattering 
channels, which are insensitive to the effects of dependent 
scattering. 

If both the atoms are located in the wavefront plane of 
the driving field, as shown in Fig. [3J the spectral depen- 
dence of the cross section is also described by two reso- 
nance features. With referring to the excitation scheme 
defined in the laboratory frame see Fig. [1] in the spe- 
cific planar geometry the double-particle self-energy part 
(|2.17[) can couple only the states |1,±1) related to ei- 
ther upper (V-type) or lower (A-type) atomic levels. As 

a consequence the resolvent operator R{E) has a block 
structure and only its 4 x 4 block, built in subspace 
|0, 0) 1 1 1 , ±1)2, |1, ±l)i|0, 0)2, can actually contribute to 
the scattering process. We subscribed the states by the 
atomic number a, b — 1,2. The eigenstates of this ma- 
trix have different parities g (even) and u (odd) reflecting 
their symmetry or antisymmetry to transposition of the 
atomic state; see Ref. [2l| . [23| . The observed resonances 
can be associated with two even-parity states symmetri- 
cally sharing the single excitation in the system of two 
atoms. Such selection rule is a consequence of the ev- 
ident configuration symmetry of the system, shown in 
inset of Fig. [3J to its rotation on any angle around the 
z axis such that the allowed transition amplitude should 
be insensitive to the atoms' positions. In contrast to the 
collinear geometry case in the planar geometry both the 
resonances have identical shapes and line widths. It is 
also interesting that for this specific excitation geometry 
both the atomic systems of either V- or A-type demon- 
strate similar spectral behavior. 

In general for random orientation of the diatomic sys- 
tem shown in Fig. 0] there are four resonances. These 
resonances can be naturally specified in the internal refer- 
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FIG. 2: (Color online) Spectral dependencies of the total cross 
section for a photon scattering from the system of two "two- 
level" V-type atoms (upper panel) and A-type atoms (lower 
panel) in the collinear excitation geometry; see inset. In the 
case of V-type atoms, in accordance with predictions of the 
Dicke model [13|, |lSj, the observed resonances demonstrate 
either super- or subradiant behavior when interatomic sepa- 
ration R becomes shorter. In the case of A-type atoms the 
resonances are less resolved and both have a line width com- 
parable with atomic natural decay rate. 



ence frame, where the quantization axis is directed along 
the internuclear axis, via following the standard defini- 
tions of diatomic molecule terms; see Ref. [2l|. There 
are two E g and terms of different parity and two dou- 
bly degenerate II g and II „ terms, which also have differ- 
ent parities. Here the defined terms are associated with 
the symmetry of the self-energy part and specified by the 
transition type in the internal frame such that the transi- 
tion dipole moment can have either projection (S term) 
or ±1 projection (II term). For random orientation all 
these resonances can be excited and in case of the V-type 
atoms the odd-parity resonances have subradiant nature 
and the even-parity ones are superradiant. In contrast in 
the case of the A-type atoms the observed resonances are 
less resolved and have comparable widths; two of them 
have rather small amplitudes (see the lower panel of Fig. 
3]). The previous configurations with the collinear and 
planar excitation geometries respectively correspond to 
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FIG. 3: (Color online) Same as in Fig. [2] but for planar 
excitation geometry. In both the excitation schemes for either 
V- or A-type atoms there is a symmetric resonance structure; 
see the text. 
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FIG. 4: (Color online) Same as in Fig. [2] but for random 
excitation geometry. For V-type atoms there are two super- 
radiant and two subradiant resonances. For A-type atoms the 
four resonances are less resolved and have line widths compa- 
rable with a single-atom natural decay rate. 



the excitations of the H g and H u , and E g and II ff , reso- 
nance pairs. Summarizing the results, we can point out 
that all the plotted dependencies demonstrate significant 
difference in the cooperative scattering dynamics result- 
ing from the similar quantum systems shown in Fig. [TJ 



B. Cooperative scattering from a collection of 
A-type atoms randomly distributed in space 

Evaluation of the resolvent operator for the situation 
of a many-particle system is a challenging task and its 
solution depends on the type of transition driven in the 
atomic system. For V-type atoms the problem can be 
solved even for a macroscopic atomic ensemble since the 
number of equations rises linearly with the number of 
atoms; see the relevant estimate given in Sec. Ill Bl In 
Ref. [18( the transformation of light scattering on macro- 
scopic atomic ensemble consisting of V-type atoms were 
analyzed as functions of the sample density. Particu- 
larly, the authors demonstrated how the smooth spectral 
dependence of the cross section, observed in the limit 
of dilute and weakly disordered distribution of atomic 



scatterers, would transform to the random speckle reso- 
nance structure in the case of strongly disordered and 
dense distribution containing more than one atom in 
the cubic wavelength. The presence of narrow sub- 
radiant Dicke-type resonance modes revealed a micro- 
cavity structure built up in an environment of randomly 
distributed atomic scatterers that can be posed as a cer- 
tain analog of Anderson-type localization of light. 

Our analysis in the previous section indicates that in 
the example of the A-type atoms the subradiant modes 
are not manifestable and such a system would be not 
suitable for observation of the localization effects. For 
coherent mechanisms of the quantum memory, which we 
keep in mind as a most interesting implication, the exis- 
tence of the localization regime would be useful but not a 
crucially important feature of the light propagation pro- 
cess. However, the spectral profile of the scattering cross 
section and its dependence on the atomic density and 
sensitivity to the level of disorder are very important, 
for example, for further consideration of an EIT-based 
memory scheme. 

Below we consider an example of the atomic system 
consisted of five A-type atoms, which is described by the 
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405 x 405 square matrix of the resolvent operator R(E). 
With evident provisoes but at least qualitatively the sys- 
tem can be considered as having many particles and can 
show a tendency toward macroscopic behavior. We show 
how the scattering process is modified when the configu- 
ration is made more dense and how this corresponds with 
the description of the problem in terms of the macro- 
scopic Maxwell theory. In macroscopic description the 
atomic system can be approximated by a homogeneous 
dielectric sphere of a small radius, which scatters light in 
accordance with the Rayleigh mechanism; see Ref. p2| . 
We fix the parameters of the dielectric sphere by the same 
density of atoms as we have in the compared microscopic 
random distribution. The calculation of the dielectric 
susceptibility were made similarly to that done earlier 
in Ref. [18j and we will publish the calculation details 
elsewhere. The key point of our numerical analysis is 
to verify the presence of the Zeeman structure, which 
manifests itself via the Raman scattering channels in the 
observed total scattering cross section. 

In Fig. [S]we show how the scattering cross section is 
modified with varying atomic density no, scaled by the 
light bar wavelength A, from no A 3 = 0.1 (dilute config- 
uration) to no A 3 = 1 (dense configuration). There are 
two reference dependencies shown in these plots and in- 
dicated by dashed and solid black curves. The dashed 
curve is the spectral profile of single-atom cross section 
(T = o" (A) multiplied by the number of atomic scatter- 
ers N = 5. The solid black curve is evaluated via the self- 
consistent macroscopic Maxwell description and repro- 
duces the scattering cross section for the Rayleigh parti- 
cle performed by a small dielectric sphere. Other depen- 
dencies subsequently show the results of microscopic cal- 
culations of the scattering cross section: (green [dashed 
light gray]) for a particular random configuration (vi- 
sualized in insets) and (red [dash-dotted dark gray]) the 
microscopic spectral profiles averaged over many random 
configurations. 

The upper panel of Fig. [5] relates to the low-density 
(i.e., dilute configuration or weak disorder) regime, which 
is insensitive to any specific location of atomic scatter- 
ers in space. Indeed the exact result evaluated with the 
microscopic model is perfectly reproducible by the sim- 
plest approximation of the cross section by the sum of 
partial contributions of all five atoms considered as in- 
dependent scatterers. This confirms the traditional vi- 
sion of light propagation through a multiparticle atomic 
ensemble as through the system of independent scatter- 
ers, which are in background of many practical scenar- 
ios of interaction of atomic systems with external fields. 
The Raman channel manifests in the scattering process, 
a direct consequence of the Zeeman degeneracy of the 
atomic ground state. In contrast, the central and bot- 
tom panels of Fig. [5] show how the scattering process is 
modified in the situation of high density and strong dis- 
order when the near-field effects are manifestable. The 
system evidently demonstrates cooperative behavior and 
the scattering mechanism becomes extremely sensitive to 
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FIG. 5: (Color online) Spectral dependencies of the total cross 
section for a photon scattering from the system of five A-type 
atoms randomly distributed in the space for different densi- 
ties n % 3 = 0.1, 0.5, 1.0. The solid gray curve indicates the 
spectral profile of the cross section for five independent scat- 
terers. The black solid curve performs the calculations for 
Rayleigh particle, which are based on macroscopic Maxwell 
theory. The green [light gray] dashed curves show the mi- 
croscopic calculations done for a particular configuration and 
the red [dark gray] dashed-dotted curves show the statistical 
averaging over many configurations. 



any specific distribution of the scatterers in space. The 
spectral profile is described by several resonances, and lo- 
cations, amplitudes, and widths are unique for each spe- 
cific configuration. However, there is a certain tendency 
to compromise the microscopically calculated scattering 
profile with the rough macroscopic prediction. The lat- 
ter keeps only the Rayleigh channel as observable in the 
self-consistent macroscopic model of the scattering pro- 
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cess. It is interesting that for any configurations, cre- 
ated randomly in the spatial distribution of the atomic 
scatterers, one of the observed resonances is preferably 
located near the vicinity of the zero detuning A ~ 0. As 
a consequence, after the configuration averaging, the sys- 
tem demonstrates scattering characteristics qualitatively 
similar to those reproduced by the macroscopic model. 



C. Application to atomic memory problem 

The considered system of A-type atoms has certain po- 
tential for light-assisted coherent redistribution of atoms 
in the ground state initiated by simultaneous action of 
strong control and weak signal modes, that is, for real- 
ization of atomic memory protocol. Let us discuss the ap- 
plicability and advantage of such a dense configuration of 
atoms for realization of light storage in atomic memories. 
At present most of the experiments and the supporting 
theoretical discussions operate with dilute configurations 
of atoms either confined with MOT at low temperature or 
existing in a warm vapor phase; see Refs. [l|-|3|- For such 
systems the standard conditions for realization of either 
EIT- or Raman-based storage schemes require an optical 
depth of around hundreds such that the macroscopic en- 
semble would typically consist of billions of atoms. The 
optimization of the memory protocol for the parameters 
of optical depth, pulse shape, etc., has been the subject 
of many discussions in literature, see Ref. [H and refer- 
ences therein. There would be an evident advantage in 
developing the memory unit with fewer atoms but with 
the same optical depth of the sample. This immediately 
readdresses the basic problem of cooperative light scat- 
tering by a dense system of the A-type configured atoms. 

The presented microscopic analysis of the scattering 
process in such systems shows that in the strong disorder 
regime the spectral profile of the cross section is generally 
described by rather complicated and randomized reso- 
nance structure contributed by both the longitudinal and 
transverse self-energy interaction parts of the resolvent 
operator. This spectrum is unique for each particular 
configuration of the atomic scatterers and has only slight 
signature of original nondisturbed atomic spectrum. This 
circumstance is a direct consequence of the complicated 
cooperative dynamics, which reflects a microcavity na- 
ture of light interaction with a strongly disordered atomic 
ensemble. 

To determine possible implications of our results to 
the problem of atomic memories, we should extend the 
presented calculations toward the ensembles consisting 
of a macroscopic number of atoms. Such an extension 
seems not so straightforward since the number of con- 
tributing equations rises exponentially with the number 
of atoms and certain simplifying approximations are evi- 
dently needed. In this sense our calculations of the scat- 
tering cross section performed for a small collection of 
atoms can be considered a precursor to calculation of the 
transmittance coefficient, which would be a key charac- 



teristic in the macroscopic description of the problem. 
Our calculations indicate preferable contribution of the 
Rayleigh mechanism in the overall cooperative scattering 
process for a density and disorder level near the Ioffc- 
Regel bound no A 3 ~ 1. It is important that in this case 
one of the absorption resonances is located in the spectral 
domain near the zero detuning for any atomic configu- 
rations and provides the desirable conditions for further 
observation of the EIT phenomenon. The presence of the 
control mode, tuned at this predictable resonance point 
and applied in any "empty" arm of the A scheme see Fig. 
[TJ would make the atomic sample transparent for a signal 
pulse. Due to controllable spectral dispersion the signal 
pulse could be delayed and effectively converted into the 
long-lived spin coherence. 

Realization of this scheme requires essentially fewer 
atoms than for dilute ensembles prepared in warm va- 
pors and in MOT experiments. Roughly for a fixed opti- 
cal depth &o ~ no% 2 L, where L is the sample length, and 
for hqX 3 -C I, the required number of atoms, allowing 
for diffraction losses, should be more than ^/^oA 3 . This 
number can be minimized if we approach the dense con- 
figuration no A 3 ~ I and make the near field effects man- 
ifestable. We are currently working on a self-consistent 
modification of the presented calculation scheme to make 
it applicable for a multiatomic ensemble and then to de- 
scribe the problem in a macroscopic limit. This can be 
done if we take into consideration the near-field effects 
only for the neighboring atoms separated by a distance of 
wavelength. For the intermediate densities with no A 3 ~ I 
we can soften our original estimate, given in Sec. Ill Bl for 
the number of equations to be solved, and can expect that 
the actual number would be scaled as d e N d™ -1 . Here 
n — 1 ~ no A performs the varying parameter denoting 
the number of the neighboring atoms, which have near 
field interfering with a selected specific atom. Our pre- 
liminary analysis shows that such a calculation algorithm 
should demonstrate a rapidly converging series with in- 
creasing n and would allow us to include the control 
mode in the entire calculation procedure. Such a mod- 
ification of the performed calculation scheme would be 
practically important and generally interesting for better 
understanding the microscopic nature of a A-type optical 
interaction in macroscopic atomic systems existing in a 
strong disorder regime. 



IV. SUMMARY 

In this paper we have studied the problem of light scat- 
tering on a collection of atoms with degenerate structure 
of the ground state, which cooperatively interact with 
the scattered light. We have discussed the difference for 
the scattering process between such system of atoms and 
well-known object of the Dicke problem, performing an 
ensemble of two-level V-type atoms. The investigation 
is specifically focused toward understanding principle as- 
pects of the scattering processes that can occur and how 
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they vary as the atomic density is varied from low values 
to levels where the mean separation between atoms is on 
the order of the radiation wavelength. For both the A- 
and V-type systems the spectral profile of the scattering 
cross section strongly depends on the particular atomic 
spatial configuration. However, in the case of the de- 
generate ground state, the presence of Raman scattering 
channels washes out visible signature of the super- and 
subradiant excitation modes in the resolvent spectrum, 
which are normally resolved in the system consisted of 
two-level atoms. We have discussed advantages of the 
considered system in the context of its possible implica- 
tions for the problem of light storage in a macroscopic 
ensemble of dense and ultracold atoms and we point out 
that the quantum memory protocol can be effectively or- 



ganized with essentially fewer atoms than in the dilute 
configuration regime. 
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